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Abstract 



An extension may be proposed to the intensity interferometer of Hanbury Brown 
and Twiss to provide the Fourier phase measurement by the use of third-order 
intensity correlations. It is well known that interferometric reconstruction of astro- 
nomical images can be obtained from second-order correlations only when a priory 
information is used about the object. The third-order intensity correlations contain 
information about the Fourier phase and need no such assumptions. In the ordinary 
way we can make measurements of the second-order intensity correlations with an 
intensity interferometer. We can also calculate the third-order intensity correlations 
from the data set of measured intensities for each distance triplet. After that we 
can reconstruct the Fourier phase from third-order correlations by using bipectral 
technique. When this is combining with the Fourier magnitude obtained from the 
second-order intensity correlations, we have the ultimate Fourier transform of the 
source brightness distribution. An inverse Fourier transform recovers the source 
image. 



According to the van Cittert-Zernike theorem any quasi-monochromatic source creates in a 
projective plane spatial coherent pattern. It is described by complex degree of coherence j(Pi , P 2 ) 
and is equal to the Fourier transform (FT) of the source brightness distribution 
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where I is the intensity of the source, A - the mean wave length, (£, 77) are coordinates on the 
source, and (p, q) are 
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Here R is the distance from the source to the points P\ and Pi in a projective plane, and {x\, yi) 
and (#2,2/2) are the coordinates of Pi and P2, respectively. Since the function 71,2 is complex, it 
has amplitude and phase. 

We can make measurements of the Fourier magnitude from the second-order intensity corre- 
lations with an intensity interferometer 
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But in so doing, information about a phase is lost. Third-order correlations provide advantages 
over second-order correlations for image recovery by allowing phase information, as reviewed by 
Lohmann and Wirnitzer (1984). 

We shall consider 1-D signal reconstruction. The extension to 2-D images is more complicated, 
but in principle, it is the same as for 1-D signal. We consider the arbitrary 1-D real signal I(x) 
whose third-order spatial correlation is given by (Webster et al. 2002, 2003) 
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Lohmann and Wirnitzer (1984) define the bispectrum B(fi, f 2 ) as the FT of I^(x\,x-2) 
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The bispectrum of I(x) may be written as 

B{h,h) = F{h)F{f 2 )F*{h + h), (5) 

where * denotes complex conjugation. We shall express the FT of I{x) in terms of its magnitude 
and phase as F(f) = A(f) ■ exp (i <f>(f)) . Analogically, the bispectrum is defined as \B(f\,f2)\ ■ 
exp (i VK/i, f 2 ))- The relationship between the Fourier phase and the bispectral phase of I{x) 
may be obtained from Eq. (5) as 
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This equation can be exploited to reconstruct the Fourier phase of I(x) from its bispectral phase. 
Writing the phases in discrete notation, allows the Fourier phase to be reconstructed, as noted 
by Bartlet et al. (1984) 
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where </>o = 0, <f>\ is arbitrary, N is the length of the sequence. Having the value of <pi, any other 
<f>k can be obtained recursively. A condition for this approaches, that <pi is arbitrary, follows from 
the fact that the bispectral phase (6) is independent of the linear phase component af. This 
gives an arbitrary offset in the reconstructed signal. 

The normalized intensity I = Al/ < I >, when used in the third-order spatial correlation in 
Eq. (2), allow the bispectral phase to be written as 

(8) 

Note that there is a sign ambiguity in the bipectral phase in Eq. (8), which results in time 
reversal indeterminacy. We can take the appropriate sign in Eq. (8) using the calculated data. 
The value of <f>\ linearly effects the phase (7); it may be often estimated from a priory data. Thus, 
by using bipectral technique, we can reconstruct the Fourier phase. When this is combining with 
the Fourier magnitude (2) obtained from the second-order intensity correlations with an intensity 
interferometer, we have the Fourier transform of the source brightness distribution (1). An inverse 
Fourier transform recovers the source image. 

Some experimental results illustrating the techniques developed above may be found in Sun- 
daramoorthy et al. (1990), Webster et al. (2002, 2003). 
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